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A proposal for a phase gate and a Mplmer-Sprensen (MS) gate in the dressed state basis is presented. In order 
to perform the multi-qubit interaction, a strong magnetic field gradient is required to couple the phonon-bus to 
the qubit states. The gate is performed using resonant microwave driving fields together with either a radio¬ 
frequency (RF) driving field, or additional detuned microwave driving fields. The gate is robust to ambient 
magnetic field fluctuations due to an applied resonant microwave driving field. Furthermore, the gate is robust 
to fluctuations in the microwave Rahi frequency and is decoupled from phonon dephasing due to a resonant RF 
or a detuned microwave driving field. This makes this new gate an attractive candidate for the implementation 
of high-fidelity microwave based multi-qubit gates. The proposal can also be realized in laser-based set-ups. 
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High-fidelity quantum gates are a crucial element in the 
growing field of quantum information processing (QIP) [1], 
Many theoretical proposals for quantum entangling gates have 
been considered for trapped ions [2-11], one of the most 
promising candidates for QIP [12-14]. These proposals have 
triggered impressive experimental realizations [15-25], Al¬ 
though these experiments within the laser based designs have 
achieved high fidelities [16, 22-25], the achieved fidelities 
within the microwave based designs have been limited. 

Considerable theoretical efforts have been made to counter 
the fidelity-damaging effects. Techniques, such as dynamical 
decoupling using echo-pulse sequences [26, 27], and its con¬ 
tinuous version using a continuously applied driving field [28- 
33], have been proposed and realized experimentally [20, 34- 
36], Recently, a combination of the continuous techniques 
with gate operators has been proposed [10, 11] and realized 
[16, 17] for the laser-induced implementations. 

In this manuscript, we introduce a scheme for a geometric 
gate (7- 0 <7-, and a M0lmer-S0rerisen (MS) gate <7 V 0 a x in 
the dressed state basis for microwave-based implementations. 
It can also be implemented in laser-based experimental set¬ 
ups. Our scheme combines the gate operator and continuous 
dynamical decoupling which results in the gate being decou¬ 
pled from the main fidelity damaging noise sources such as 
ambient magnetic field and Rabi frequency fluctuations. 

In what follows, we derive the gate Hamiltonian for a 
microwave-based design, while highlighting the dynamical 
decoupling processes which protect from the main noise 
sources. We then go on showing explicitly how our scheme 
could also be implemented using laser rather than microwave 
radiation. 


I. METHOD 

In the microwave-based implementation the required two 
spin states defined as {11), |0)} are separated by an energy 
splitting in the microwave regime [20, 21, 34]. The spin-spin 
interaction, which constructs multi-qubit gates, is induced via 
the exchange of virtual phonons between the ions, in the weak 
coupling regime [3], or real phonons in the strong coupling 
regime [4, 8]. Therefore, a coupling between the spin and the 


interaction-mediating phonon is needed. Due to the negligi¬ 
bly small Lamb-Dicke parameter associated with microwave 
radiation, inducing a sufficiently strong spin-phonon coupling 
requires the use of a large magnetic field gradient. 

There are two basic approaches to obtain a magnetic field 
gradient. One approach makes use of an oscillating magnetic 
field gradient induced by the microwave driving field in the 
near field regime, where the spin-phonon coupling is already 
linked to the driving fields which perform the gate [9, 21], 
This enables to utilize qubits that can have first order mag¬ 
netic field insensitive transitions [22] ; namely, the qubit states 
should be magnetic field sensitive, yet the energy difference 
between the two states is insensitive to the first order in the 
magnetic field. Therefore, these qubits are decoupled from 
the ambient magnetic noise, and thus can be considered as 
clock qubits. 

The other approach uses a static magnetic field gradient 
[37, 38], Here, the qubits must have first order field sensi¬ 
tive transitions, making them extremely sensitive to the ambi¬ 
ent magnetic noise, which imposes the use of dynamical de¬ 
coupling techniques. Furthermore, in this approach the spin- 
phonon coupling is not linked to driving field transitions that 
eventually perform the gate. In order to analyze this setting, 
one can apply a polaron-like transformation [39]. In that case, 
there always remain undesired transitions, which originate di¬ 
rectly from the static magnetic field gradient [40, 41], Not 
accounting for these transitions may damage the gate fidelity, 
and adiabatically eliminating them imposes restrictions on ex¬ 
perimental parameters, such as an upper bound to the effec¬ 
tive Lamb-Dicke parameter, which can significantly increase 
the gate duration. Instead of using the dressing fields as the 
gate performing term, in this proposal, the gate is generated 
directly by the static magnetic field gradient. 

To realize a gate in the strong coupling regime, it is nec¬ 
essary to have sideband addressing of the resolvable secu¬ 
lar frequencies, namely to couple to only one common mode 
of motion of secular frequency v. This is done by applying 
a microwave driving field on-resonance with the bare qubit 
{11), |0)}, whose Rabi frequency is nearly equal to the se¬ 
lected common mode of motion, D = v — e, and far detuned 
from all the other modes. For a two-qubit gate, an increase of 
the gate speed could be obtained by coupling to two modes. 
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with opposite detunings. This is done by setting the mi¬ 
crowave Rabi frequency exactly between the two modes; how¬ 
ever, for the simplicity of the derivation, we only consider 
coupling to one mode. 


II. DERIVATION OF THE HAMILTONIAN 

Due to the constant magnetic gradient and the difference in 
the magnetic moments of the spin states |0) and |1), a state- 
dependent force is obtained. In addition, each ion feels dif¬ 
ferent magnetic field, and as consequence, different Zeeman 
splitting. Therefore, the Hamiltonian can be written as 

H = vb"b + ^ ~ 2 ~ a z + ~ 2 ~ + °z + ^ a x cos > 

(1) 

which includes the selected vibrational mode with secular fre¬ 
quency v, the bare state energy structure of the j th ion, the 
magnetic field gradient, and resonant microwave driving fields 
applied to the j th ion respectively (fig. 1). The Pauli spin ma¬ 
trices <7- = 11} (11 — |0) (0|, CJ T = |1) (0| + |0) (11, and b are 
the phonon creation and annihilation operators, respectively, 
and r\ J = gji « d z BZJ / \/2m V 3 i s the effective Lamb-Dicke pa¬ 
rameter, with ZJ being the selected normal mode coefficient. 
For simplicity we assume that ZJ is the same for all the ions, 
namely we are coupled to the center of mass mode. 

Large magnetic gradient enables single addressing [42], 
where a microwave driving field that is on resonance with one 
ion transition, is far detuned from the other ion transitions. 
The neglected off-resonant driving fields contribute an A.C 
Stark shift which results in an undesirable phase shift, and 
will be discussed in the next section. 

According to the microwave driving field’s initial phase, 
the dressed states are determined. For simplicity, in the in¬ 
teraction picture we consider a vanishing initial phase, such 
that the dressed state basis is the eigenstates of o x , namely 

l M ) = (|1) + |0)) /a/2. \d) = (|1) — |0)) /\/2 (1). 

Moving to the interaction picture with respect to the mi¬ 
crowave energy splitting, £,-(c0q/2)c 7/, followed by trans¬ 
forming to the dressed state basis, the operators are trans¬ 
formed as follows: a v —> S z ,O y —> S y , and <7- -A —S x , by where 
S a are the Pauli matrices in the a direction, in the dressed state 
basis. Thus, the Hamiltonian of eq.l is transformed and can 
be written as 

H I = vtfb + Y.[- V ~Y^+b)Si+^Si S j, (2) 

after making the rotating wave approximation (RWA), when 
assuming £2 <C 4 C 0 q. The neglected fast rotating terms of the 
microwvae driving fields contribute an A.C Stark shift which 
results in an undesirable phase shift, and will be discussed in 
the next section. 

In the interaction picture with respect to the dressed state 
energy splitting, (£2/2 )Y,jSl, the red sideband transition. 
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Figure 1. Moving from the bare states to the dressed states. In 

the microwave-based designs, due to the magnetic gradient, each 
ion feels different magnetic field, and therefore, it has a different 

bare state energy splitting J_j ai. By applying microwave driv¬ 
ing fields on resonance with the bare state energy splitting of the 
different ions, we move to the dressed state basis, in a perpendicu¬ 
lar direction y . The energy gap in the perpendicular direction 
protects against the magnetic noise. 

In the laser-based designs, if the qubit transition frequency is in the 
optical regime, instead of the resonant microwave driving fields, we 
may apply a single resonant laser transition, in order to move to the 
dressed states. 


which originates directly from the magnetic field gradient, can 
be written as 


Hred = - ^ £ ( fct e ‘ Vl +h.c) ( s j + e int +h.c). (3) 

In the second-order perturbation approach the red sideband 
transitions result in a flip-flop Hamiltonian Y,i,jS ] + S’_ + h.c 
with an undesired AC Stark shift term, which is coupled to the 
number of phonons £ ; S{b'b, behaving as a phonon dephasing 
source. Instead of adiabatically eliminating these transitions, 
we use them as our gate-performing terms, while suppressing 
the dephasing term. Transforming eq. 3 to the desired gate 
transitions can be achieved in two ways: (a) by applying a sin¬ 
gle RF driving field (fig. 2), which is on-resonance with the 
dressed state energy splitting determined by the microwave 
driving field, and polarized in the axial direction, which is de¬ 
termined by the external magnetic field, (b) By applying an¬ 
other microwave driving fields (fig. 3), which are detuned 
from the bare state energy splitting by exactly the Rabi fre¬ 
quency of the resonant microwave driving field, 8=0., and 
phase-locked with it. These additional driving fields can be 
represented as 


( a)H r f = £2, £ ai cos £2r 

. 1 (4) 

=0 2 J^ &x cos ((Oq - £2 )t 
j 

Following the same stages of derivation up to this point, 
namely, moving to the interaction picture with respect to the 
qubit splitting, followed by moving to the interaction picture 
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Figure 2. Moving from the dressed states to the double dressed 
states using RF held. By applying a single RF driving field on res¬ 
onance with the dressed state energy splitting we move to 

the double-dressed state basis, in a perpendicular direction % Y.j ■ 
The energy gap in the perpendicular direction protects against the 
noise that originates from the Rabi frequency fluctuations £2. 


according to the dressed state energy structure, the additional 
terms (eq.4a,b) can be written as 


o ' (5) 

j 

after making the RWA, where we assume that £2, ,£2o <C 4£2. 
The neglected fast rotating terms contribute an A.C Stark shift 
which results in an undesirable phase shift. Using detuned mi¬ 
crowave fields instead of a single RF field, contribute another 
A.C Stark shift, due to the neglected off-resonant contribution. 
These undesired A.C Stark shifts will be discussed in the next 
section. 

Transforming to the double-dressed state basis, and moving 
to the interaction picture with respect to the time-independent 
term of the RF driving field or the detuned microwave driving 
field (eq.5), the red sideband transitions (eq.3) become 


Figure 3. Moving from the dressed states to the double dressed 
states using detuned transition in the bare state basis. In the 

microwave-based designs, due to the magnetic gradient, each ion 
feels different magnetic field, and therefore, it has a different bare 

state energy splitting ^ L / O') ■ Applying another microwave driv¬ 
ing fields, which are detuned from the bare state energy splitting by 
exactly the Rabi frequency of the resonant microwave driving field, 
8 = £2, and phase-locked with it, is effectively equivalent to driv¬ 
ing resonant transition with the dressed state energy structure. This 
takes us to the double-dressed state basis, in a perpendicular direc¬ 
tion D ~/ J_j a/. The energy gap in the perpendicular direction protects 
against the noise that originates from the Rabi frequency fluctuations 
£2. A single detuned microwave driving field can also be used for 
the laser-based designs, if the qubit transition frequency is in the mi¬ 
crowave regime. If the qubit transition frequency is in the optical 
regime, instead of the detuned microwave driving field, we may ap¬ 
ply a single laser transition, which is detuned from the bare state en¬ 
ergy splitting by exactly the Rabi frequency of the resonant laser tran¬ 
sition, 8 = £2. In that way the system moves to the double-dressed 
states which are decoupled from the Rabi frequency fluctuations of 
the resonant laser transition. 


eq.7. The Magnus expansion can then be written as 


1 

2 


U{t) =exp(—/ [ dt'H gate (t') 
Jo 

fdt'f dt" [H gate (t'),H gate (t")\ + ...). 


( 8 ) 


H'j = ^ (bV v ' + h.c ) ^ (cs ] z - <xf+ a!_e ,art ^j e ,nt + h.c 

( 6 ) 

where for simplicity we have merged the notations: £2, = 

£ 12 / 2 . 

By setting e := v £2 £2,, the only transition which is 

not suppressed, and thus survives, is the desired phase gate 
transition in the double-dressed state basis, or the MS gate 
transition in the dressed state basis 


'For this specific Hamiltonian, only the first two orders of the 
Magnus expansion contribute, whereas higher orders vanish 
[44], An alternative approach to obtain the same result is to 
follow the MS gate notation in Ref. [4], According to it, the 
gate Hamiltonian can be represented by 

77 V 

Hgate = {xcoset + psinet) =: f(t)x + g{t)p. (9) 

v8 j 


H S are =7^ {b'eP+be-**) = Y (b'eP+be-**) 

(7) 

In order to obtain the effective Hamiltonian of the gate, we 
can use the Magnus expansion [43] of the time propagator 
corresponding to the time-dependent Hamiltonian shown in 


where x = {b^ +b)/y/2 and p = —i(b ^ —b)/\fl are the dimen¬ 
sionless operators of the position and momentum respectively. 
Since the first two terms of eq.9 commute with their mutual 
commutator, we obtain the gate’s exact time propagator given 
by 

U[t)=e~ iA{t) e- iF(t)x e- iG{ ' t)p , 


( 10 ) 
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where 

F (t) = f(t')dt' = -j= o'sine* 


We have dropped the fast rotating terms of the second driv¬ 
ing field (the RF field or the additional detuned microwave 
(11) driving fields). Therefore, in eq. 5, for the RF driving case, 
we should have additional rotating terms: 


GO) = ^ g(t r )dt' = Y IX J 0 “ cos£t ) G2) 

Mt) = ~j 0 = -~^Y i O l z al ^-^sin2e^ 

(13) 

In order to obtain the required behavior of the gate, we have 
two restrictions: First, we should suppress the entanglement 
of the spin to the motion, that is to say, we set F(x) = G(z) = 

0, by choosing T = 2 itK/e, with an integer K, representing 
the integer number of circles in phase space. Second, we 
would like to have a maximally entangled state when start¬ 
ing with a separable state, so we require A{x) = —(k/%)gIg1. 
These constraints determine the connection between the de¬ 
tuning and the Rabi frequency of the sideband transitions, and 
the gate time, £ = !] v\fK and T = 2K\/~K/r ]\\ respectively. 

All of the assumptions made during the above derivation 
can be summarized as 

i]V /4\/K = e/4 <C O r /4 <C v ~ £2 <C 4<a 0 . (14) 

in. CALCULATING THE NEGLECTED TERMS 

As was mentioned above, the neglected off-resonant mi¬ 
crowave driving fields in eq. 1, and the fast rotating terms 
of the second driving field (the RF field or the additional de¬ 
tuned microwave driving fields) in eq. 5 contribute an A.C 
Stark shift which results in an undesirable phase shift. In ad¬ 
dition, except from the gate operator eq. 7, we have also ne¬ 
glected other terms from eq. 6. In this section we will estimate 
the leading contribution of these fidelity damaging terms, and 
suggest a way to counter them using spin-echo. 

large magnetic gradient enables single addressing [42], 
since each ion feels different magnetic field, and as conse¬ 
quence, a different Zeeman splitting. Suppose we have two 
ions in our trap, aligned in the z direction. Thus, the energy 
difference between the two ions is 

Amo = 0J,| ft),] = gfl b d z BAZ, (15) 

where 

/ 2e 2 \ 1/3 

AZ = -—- (16) 

\AtzeqMv- J 

is the distance between the ions in the axial direction, £o is the 
electric permittivity, M is the ion mass, and e is the elementary 
charge. Taking the off-resonant microwave driving fields into 
account, another term appears in the rotating frame of the bare 
energy structure 

j (ole**# + G 2 e ~ ih ^ + h.c. (17) 


- Y( s ++ s l) e 2 iar + h -c- ( 18 ) 

In addition, in the previous section we have also neglected 
the following terms from eq. 6: 

^ (b f e iv, +h.c) ^(^-G J + e i£lrt + G J _e~ iartS je int , 

(19) 

and 

a/(*V (v+0)t + A.c). (20) 

We can evaluate the leading contributions of these terms by 
taking into account the time-independent terms of Magnus ex¬ 
pansion. For eq. 17,18, the leading contribution comes from 
the second order of Magnus expansion, and gives the follow¬ 
ing A.C Stark shift: 


~O r 


^) 2 


O r 

40. 



The leading contributions of eq. 19 are: 


2 

Or 



( 21 ) 


( 22 ) 


which is an A.C Stark shift also coupled to phonons, and 


- 


1 


V 4 J \0 2 v + £2 


(o]G 2 +h.c), (23) 


which couples between the two ions. Another coupling term 
comes from the leading contribution of eq. 20: 


-2 (m : 


2 1 ,2 

-<7 (7 

v 4 ) v+n z z 


(24) 


The A.C Stark shifts of eq. 21 and 22 cause dephasing, and 
damage the gate fidelity. In order to compensate for the de¬ 
phasing A.C Stark shifts, we can move to the interaction pic¬ 
ture with respect to these terms, by applying a global rotation 
about the z axis. Alternatively, this undesired contribution can 
be compensated for by applying a single 7T-pulse after half of 
the gate time. Note that only when an integer number of cy¬ 
cles in phase space have been completed, should a 7T-pulse 
be taken; therefore, a 7 T-pulse in the middle of the gate deter¬ 
mines the total number of cycles, K, to be even. 

Instead of applying a 7T-pulse, we can use an alternative 
procedure, by applying a 7T-phase flip of the RF driving field, 
in the middle of the gate [8, 45]. This is due to the following 
three reasons: 1.) Applying a /r-phasc flip of the RF driving 
field, does not change the gate Hamiltonian (eq. 7). 2.) It turns 
out that all the leading contributions of the A.C Stark shifts 
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(eq. 21,22) are an odd function of £l r , therefore, changing 
the sign of £2 r , would change the sign of the undesired terms. 
3.) These A.C Stark shifts also commute with the gate Hamil¬ 
tonian (eq. 7), and therefore, one dynamical-decoupling-like 
7T-phase flip, applied in the middle of the gate, is sufficient to 
counter the undesired A.C Stark shifts. Another advantage of 
applying a 7T-phase flip of the RF driving field is that it can 
work even for K = 1, as opposed to the 7T-pulse technique. 

Since eq. 22 is also coupled to the phononic states /?7?, 
which are populated during the gate duration, (especially due 
to the electric field effect, as will be discussed below), apply¬ 
ing a single n —phase flip may not be sufficient to obtain high 
fidelity. For that reason, we can apply a sequence of K —phase 
flips. Note that this sequence can be performed even within a 
single cycle in phase plane (K = 1). 

We can evaluate the higher order contribution of the ne¬ 
glected terms using the simulation of the gate performance, 
including all the terms mentioned above without any approx¬ 
imation. In the simulation we have two Yb ions of mass 
m = \13amu, magnetic gradient d z B = 65 T / m , secular fre¬ 
quency v = 2 K ■ 500 kHz, effective Lamb-Dicke parameter 
T] = 0.01, microwave field Rabi frequency £2 = 2k ■ 495 kHz, 
RF field Rabi frequency £1, = 2k ■ 99kHz, one loop in phase 
space k = 1, the off-resonant microwave driving field con¬ 
tribution where Afflo = 2 n ■ 5 MHz, starting in the vibrational 
ground state, namely, the initial number of phonons is N = 0. 
This gives a gate duration of t gate = 0.2 msec. In addition, we 
applied the /r-phase flip in the middle of the gate (fig. 4). 

We find that the higher order expansion of eq. 17, is an 
even function of £l r , which is not refocused by the spin-echo. 
It yields infidelity of IF = 0.002 (fig. 5) . However, although 
the current set-up does not allow it, in future implementations, 
the ions may be aligned along a perpendicular direction to 
the magnetic gradient, such that the ions would feel the same 
magnetic field, and thus have the same bare energy structure 
due to Zeeman splitting. Therefore, the undesired contribution 
of eq. 17 would vanish. Simulating all the terms but eq. 17 
yields infidelity of IF = 3 • 10 -4 (fig. 5), getting closer to the 
fault tolerance threshold. 

For these suggested parameters, the contributions of the 
coupling terms, eq. 23 and 24, correspond to infidelity of 
1.2 • 10~ 5 and 5 • 10~ 5 respectively. 


IV. ROBUSTNESS TO NOISE 

There are several noise sources in this proposal, which are 
significantly suppressed using concatenated continuous dy¬ 
namical decoupling [31], The main noise source originates 
from ambient magnetic field fluctuations. This is suppressed 
by the resonant microwave driving field [20, 34], which can 
be readily shown in the dressed state basis, the eigenstates of 
a x . The x direction is orthogonal to the quantization axis, and 
therefore to the magnetic noise which is described by the op¬ 
erator <7- (fig. 1); thus, by moving to the interaction picture 
of the dressed state energy structure, as was done in eq.3, the 
magnetic noise is suppressed in the RWA, exactly as is the 
case in the gate scheme proposed by Bermudez et al. [10]. 
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Figure 4. Entangling gate performance. The dynamics of two qubit 
gate starting with \dd) as the initial state with the ground state of the 
vibration mode is plotted. We obtain the maximum entangled state 
IT-i-) = (| dd) + i\uu))/ \/2 at the end of the gate, after applying a K- 
phase flip of the RF driving field in the middle of the gate. Counting 
from above at time = 0.08 msec, the first (cyan) curve is P\dd)Uld\- the 
second (blue) is the fidelity of being in the maximum entangled state 
F = {*P+|p | l F+), the third (red) is the imaginary part of P\dd)(uu\’ 
the fourth (cyan) is P\ uu \/ uu \, and the lowest one (green) is the real 
part ofp| rfd) { 

uu \* 



time[o.ims] 


Figure 5. Zooming-in of the gate fidelity (eq. 4) assuming differ¬ 
ent terms in the Hamiltonian. From the comparison between the 
achieved fidelity of the case where the electric field and the single ad¬ 
dressing vanish (pink asterisks), to the case where only the single ad¬ 
dressing is taken into consideration (cyan circles), we understand that 
the infidelity due to the single addressing is IF = 0.002. Without any 
assumption, while applying only a single 7T~phase flip in the middle 
of the gate (blue triangles), the gate fidelity is F = 0.995. By increas¬ 
ing the number of n —phase flips the fidelity is increased as-well. In 
the case where we apply 19 K —phase flips during the gate (red rect¬ 
angular) the achieved fidelity almost reaches F = 0.997. When we 
suppress the single addressing (in future configurations) and apply 
99 K —phase flips (green pluses) the infidelity is IF = 3e~ 4 , almost 
reaching the fault tolerance threshold. 


The second important noise source originates from fluctua¬ 
tions in the Rabi frequency of the resonant microwave driving 
field. The Rabi frequency fluctuations are one of the most 
fidelity damaging sources in reference [16]. To protect our 
gate from this, we use a resonant RF or a detuned microwave 
driving field (eq.5), whose effect can be shown in the double- 
dressed state basis, the eigenstates of a z . The z direction is 
orthogonal to the resonant microwave driving field’s polariza¬ 
tion, and therefore to the Rabi frequency fluctuations, which 
is described by the operator cr v (fig. 2, 3); thus, by moving 
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to the interaction picture of the double-dressed state energy 
structure, as was done in eq.7, the suppression of the Rabi fre¬ 
quency noise is achieved. Since the Rabi frequency of the res¬ 
onant microwave field is much larger than the Rabi frequency 
of the RF or the detuned microwave driving field, the noise 
originated from fluctuations in one of those driving fields is 
less fidelity-damaging. 

We can estimate the infidelity due to these noise sources. 
Previous experiments [20, 34] have shown that the dephasing 
rate of the ambient magnetic field noise was suppressed by 
the dressing microwave radiation to less than Snis(2()kHz) = 1 
Hz. With higher dressing Rabi frequencies the dephasing rate 
would be decreased as 


S B j(n) = S B . B (20kHz) =0.0016 Hz, (25) 

assuming that the correlation time of the magnetic noise is 
much longer than 50 flsec. 

The dephasing rate due to the Rabi frequency fluctuations 
without additional dressing fields can be approximated as 
.S'o o(O) k, 0.010. The time correlation of this noise source 
is long in most experiments, i.e., of the same order as the two 
qubit gate rate. Thus, by applying the second dressing fields, 
the first order noise is completely suppressed Sq o(Qi ) ss 
0. Yet, the second order contribution of this noise, which 
is an A.C Stark shift, becomes dominant and behaves as 
Sq ci(0) 2 /£!, . The second order contribution of the Rabi fre¬ 
quency noise source has a long time correlation, and since it 
is an odd function of Q. ri with the spin-echo n —phase flip, its 
contribution is negligible. 

The dephasing rate due to the Rabi frequency of the last 
dressing field, the RF or the detuned microwave driving field, 
can be approximated as before as Sa r ,o. r {0) ~ 0.01 H r . We 
can assume that the time correlation of this noise source is 
very long, therefore, it can be refocused using the spin-echo. 

In total, the infidelity caused by the ambient magnetic field 
and the fluctuations of the first and second dressing fields can 
be estimated as IF = 3 • 10 7 for the 0.2 msec gate duration. 
This is much smaller than the fault tolerance threshold, and 
therefore the gate is robust to the main noise sources. 

High fidelity depends on good timing (fig. 5). In this re¬ 
gard, techniques such as pulse shaping as were suggested by 
C. F. Roos [8], and Hayes et al. [45] have been developed in 
order to overcome this problem. Yet, these techniques can not 
be implemented in the proposed scheme for the microwave- 
based implementations, since we can not modify the ampli¬ 
tude of the gate operator (vrj/4) (eq. 7), as these techniques 
require. Nevertheless, microwave signals can be created us¬ 
ing highly accurate arbitrary waveform generators which can 
address this timing issue. 

In addition, timing depends also on the stability of the Rabi 
frequency of the resonant microwave H and secular frequency 
V (eq. 13). The fluctuations in the Rabi frequency are sup¬ 
pressed by the second dressing fields, and as is estimated 
at [22], the motional error of should be less then 10 3 in a 
0.2 msec gate duration. 


V. UNDESIRED EFFECT OF THE RF DRIVING FIELD 


Applying a driving field closely-detuned from the selected 
mode of motion coherently increases the number of phonons 
in the system. Although the driving field is magnetically cou¬ 
pled to the spin, its associated electric field is coupled directly 
to the motion, and causes displacement of the ion motion. 
This interaction can be written as 


H c i = eE -z = £Ie {b^ +b) cos(v — e)t, 


(26) 


where FI/, = 


eE z 

\Jlmv 


and e is the ion charge. Fortunately, the de¬ 


tuning e reduces these effects, and protects against excessive 
motional displacement. In the motion’s interaction picture, 
the displacement becomes H^, = ( b T e' £t +h.c ), after the 
RWA, when assuming He <C 4v. Since the gate Hamiltonian 
and the additional undesired displacement term have the same 
rotating exponents, the Hamiltonian of both of them can be 
represented by 


»-.= (7E”! + t) + (27) 

where, as before, an exact time propagator exists. Complet¬ 
ing an integer number of circles in phase space results in the 
desired phase gate transition with a correction which can be 
treated by moving to the interaction picture, or by either a 
7T-pulse or a /r-phase flip after half of the gate time, as was 
discussed above. 

The real influence of the electric coupling is the number of 
phonons that are coherently added during the gate operation, 
before closing a circle in phase space. It can be easily un¬ 
derstood in the representation of the Magnus expansion. The 
first order of the Magnus expansion (eq. 8) is the displace¬ 
ment operator, which is maximized in half of a cycle in phase 
space. For example, at t = z/2K, the displacement opera¬ 
tor is D(a( 2 ^)) = e“(x) fc+ - ,! - c , where cc(^) = 

If we initialize our system in the motional ground state, af¬ 
ter t = t/2 K, the motional state becomes the coherent state 
\y)^ = |ct( 2 ^ 0 ), whose average number of phonons could 

reach < N >= |a(^)|“ ~ (He/e) 2 , if e -c He- If there are 
too many phonons in the system important assumptions such 
as the Lamb-Dicke assumption <v!] \/N /4 <C H,±t:) made in 
eq. 6 would not hold and lead to a reduction in the achievable 
fidelity. 

It is therefore advantageous to ensure that He "C £, such 
that the electric coupling can be adiabatically eliminated in the 
second-order perturbation approach. To fulfil the restrictions 
in eq. 14 (e <C H r ), the ratio between the magnetic coupling 
and the electric coupling must be large. If it is not the case, 
and the electric field is large, a sequence of many k —phase 
flips can be applied in order to suppress the population of high 
phonon states. Note that this sequence can be performed even 
within a single cycle in phase plane (K = 1). 

In the near field regime, the magnetic field dominates 
over the electric field by a factor of the wavelength divided 
by the size of the system. Therefore, the ratio between 
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the magnetic coupling and the electric coupling can be es¬ 
timated in the following way. For a RF field of frequency 
citleastoneofthespinstatesshouldbeKHz and dimension 1cm, 
the typical velocity of the electrons is 3 x 1 0 (, c. where c is 
the speed of light. In that case the electric term in the Hamil¬ 
tonian would be £Ie =evxBxx o, where e is the ion charge, B 
the magnetic field and xq is the standard deviation of the wave 
function which is of the order on 10 nm. Thus, the ratio be¬ 
tween the magnetic and the electric energy is ^ ss 30. 

By reducing the dimension of the system a better reduction 
factor can be achieved. It means that even in the worst case, 
when lacking the ability to set the polarization of the RF driv¬ 
ing field along the axial axis, the gate could still be performed 
with high fidelity. 

In addition to the above, unwanted E-fields can be created 
due to the finite resistance of a coil which creates the required 
RF field. In order to investigate this we simulated a coil that 
consists of 3 turns of copper wire with a wire diameter of 6.7 
mm and a coil diameter of 4 cm. The coil is placed 1 cm from 
the ions. To make this even more accurate we have included 
the two feed-lines going to the coil. The E-field along the trap 
axis created due to the potential drop across the coils oper¬ 
ated at a frequency of 500 kHz and a current sufficient to ob¬ 
tain LI, = 2k ■ 100 kHz has been found to be 2.04 • 10~ 5 V/m, 
which is negligible compared to the Lorentz force discussed 
above. 

Adding the rotating electric field into the simulation of the 
gate performance, with £Ie = Ll r / 30, yields infidelity of IF = 
5 ■ 1 0 :i after a single spin-echo (fig. 5). However, assuming 
that the two ions have the same bare energy gap (such that eq. 
17 can be dropped), after a sequence of 19 and 99 phase flips, 
the infidelity reduces to 3 • 10 3 and 3 • 10 4 respectively (fig. 
5). 


VI. MOVING TO THE INTERACTION PICTURES 

In our proposal we move three times to different inter¬ 
action pictures, with respect to: (1) the bare energy struc¬ 
ture Y,j( ® o / 2 )< t /, ( 2 ) the dressed state energy structure 
(£1/2) Y., <? A / , and (3) the double dressed state energy struc¬ 
ture (£1,-/2) Y./ o/ . The way it is done in the lab is the follow¬ 
ing. The driving fields accumulate the phases corresponding 
to the first interaction picture with respect to the bare energy 
structure, namely, the system automatically moves to the first 
interaction picture. 

As regards to moving to the second interaction picture with 
respect to the dressed state energy structure, at the end of the 
experiment (after time t), after eliminating the RF field, we 
operate only with the microwave dressing fields for additional 
time t a dd, during which the accumulated phase of the second 
interaction picture trivially vanishes, namely £1(t + t a dd )/2 = 
2k n for an integer n. However, since we cannot shut down 
the static magnetic field gradient, we have to change the mi¬ 
crowave Rabi frequency (fi —>• £l new ), such that it would be far 
detuned from the secular frequencies, thus adiabatically elim¬ 
inating the magnetic gradient’s contribution. Now the condi¬ 


tion should be (£1/ 2)t+ (£l new /2)t a dd = 2 Kn for an integer 
n. 

The third interaction picture with respect to the double 
dressed energy structure, is carried via the K— phase flip of 
the RF driving field. 

VII. REALIZATION WITH LASERS 

Until now, we have discussed the realization of the gate in 
a microwave-based implementation. However, with a similar 
derivation, our scheme could also be implemented in laser- 
based set-ups. Dealing with lasers, the spin to phonon cou¬ 
pling is carried via the laser’s sufficient large momentum, k. 
Therefore, the ion’s two spin states no longer have to be sen¬ 
sitive to the magnetic field. Here, we may use trapped ions, 
with a qubit transition frequency in the optical, or even in the 
microwave regime, using resonant or Raman transitions, re¬ 
spectively. 

After the RWA, where we assume that the Rabi frequency 
of the resonant laser transition (fig. 1), or the effective Rabi 
frequency of the counter-propagating Raman beams (fig. 6), 
is much lower than the qubit energy splitting, £1 <C 4oJo, the 
Hamiltonian of the ion motion and the laser interaction in the 
rotating frame of the bare state energy structure (c %/2)o- is 
given by 

Hj = vb'b + j Ya {^{e ,riL ^ +b) + h.cj , (28) 

where /]/ = kZ/\/2mv is the Lamb-Dicke parameter, with Z 
being the normal mode coefficient. 

Performing the experiment in the Lamb-Dicke regime, 
where v\lVN <C 1, and N is the number of phonons in the sys¬ 
tem, enables one to expand the Lamb-Dicke exponent to the 
second order in )]/ . Thus, we obtain a similar Hamiltonian 
to the Hamiltonian of the microwave-based implementation 
(eq.l) given by 

Hi = vb' b 

+ J L { 1a x + + b ) ~ Y (* f + b) (b t + b) |, 

(29) 

where the carrier transition performs the dynamical decou¬ 
pling, the sideband terms would perform the gate, and the ad¬ 
ditional term is an undesired phonon dephasing term. 

After transforming to the dressed state basis, and moving 
to the interaction picture with respect to the dressed state en¬ 
ergy splitting, the sideband transitions, similar to those from 
eq.3 in the equivalent microwave case, are observed with the 
undesired phonon-dephasing term [46, 47]. This is given by 

H,ed = ~i -^ —/t.c'j ( b r e ,vt +h.c ) 

2 ' ' ' ( 30 ) 

— “^4 (b^e m + h.c) (b'e' vt +h.c). 

j 
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Figure 6. Moving from the bare states to the dressed states using 
Raman beams. In the laser-based designs, if the qubit transition fre¬ 
quency is in the microwave regime, hy applying counter-propagating 
Raman beams that effectively generate a resonant transition of the 
bare state basis, we move to the dressed state basis, in a perpen¬ 
dicular direction ^ Y.j ■ We demand the counter-propagation of 
the beams in order to obtain the coupling to the phonons, which is 
needed for the entangling gate. 


Transforming eq. 30 to the desired gate transitions, while 
suppressing the dephasing term, can be achieved in several 
ways. One could apply a resonant RF driving field (fig. 2), 
as was done in the microwave based design discussed above. 
If the bare state energy splitting is in the microwave regime 
one can apply an additional detuned microwave driving field 
(fig. 3), which can alternatively be realized by applying co¬ 
propagating Raman lasers to perform the detuned transition as 
shown in fig. 7. If the bare state energy splitting is in the opti¬ 
cal regime, we could use a detuned laser driving field, whose 
additional spin-phonon coupling terms are adiabatically sup¬ 
pressed (fig. 3). For simplicity, we consider the case of ap¬ 
plying an RF driving field given by H r f = Q-rlLj a z cos Or, as 
was done in the microwave based scheme described above. 

In the interaction picture according to the dressed state en¬ 
ergy structure, the RF driving field becomes (— £l r /2)Y,jS 1 x , 
as shown in eq.5(a). Transforming to the double-dressed state 
basis, and moving to the interaction picture with respect to 
the new energy structure, the noise, which comes from fluc¬ 
tuations in the Rabi frequency induced by the laser driving 
fields, is suppressed and we are left with the phase gate in the 
double-dressed state basis, or the MS gate in the dressed state 
basis which is given by 

H gate = £ a!(lSe iEt +h.c) = e iet +h.c) 

(31) 

respectively and is equivalent to the microwave-based case. 

VIII. SUMMARY 

We have presented a multi-qubit gate proposal for trapped 
ions which can be used in both, microwave-based and laser- 
based, implementations. The gate is robust to the main fi- 


Figure 7. Moving from the dressed states to the double dressed 
states using co-propagating Raman beams to generate a detuned 
transition in the bare state basis. In the laser-based designs, if the 
qubit transition frequency is in the microwave regime, applying an¬ 
other co-propagating Raman beams that generate the detuned transi¬ 
tion in the base state basis, such that 8 = £2, is effectively equivalent 
to driving resonant transition with the dressed state energy structure. 
This takes us to the double-dressed state basis, in a perpendicular 
direction ai . The energy gap in the perpendicular direction 

protects against the noise that originates from the Rabi frequency 
fluctuations £2. 


delity damaging noise sources (ambient magnetic field fluctu¬ 
ations and Rabi frequency fluctuations) and is decoupled from 
phonons. Our proposal builds on the scheme proposed by 
Lemmer el al. [11], however, our gate can be implemented 
using microwave as well as laser based experimental setups. 
The main difference is that in our proposal the gate transitions 
originate from the resonant driving field, such that we can use 
the carrier transition for performing continuous dynamical de¬ 
coupling, instead of generating it with an additional driving 
field. We use a second driving field for both, transforming 
the red-sideband transitions to the desired gate transitions and 
for continuously applying additional dynamical decoupling to 
suppress noise from Rabi frequency fluctuations, which orig¬ 
inate from the first resonant driving field. 
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